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INTRODUCTION 


The Institute for Computer Applications in , Science and Engineering 
(ICASE) is operated at the Langley Research Center (LaRC) of NASA by the 
Universities Space Research Association (USRA) under a contract with the 
Center. USRA is a nonprofit consortium of major U. S. colleges and 
universities. 

The Institute conducts unclassified basic research in applied 
mathematics, numerical analysis, and computer science in order to extend and 
improve problem-solving capabilities in science and engineering, particularly 
in aeronautics and space. 

ICASE has a small permanent staff. Research is conducted primarily by 
visiting scientists from universities and from industry, who have resident 
appointments for limited periods of time, and by consultants. Members of 
NASA's research staff also may be residents at ICASE for limited periods. 

The major categories of the current ICASE research program are: 

a. Numerical methods, with particular emphasis on the development and 
analysis of basic numerical algorithms; 

b. Control and parameter identification problems, with emphasis on 
effective numerical methods; 

c. Computational problems in engineering and the physical sciences, 
particularly fluid dynamics, acoustics, and structural analysis; 

d. Computer systems and software, especially vector and parallel 
computers . 

ICASE reports are considered to be primarily preprints of manuscripts 
that have been submitted to appropriate research journals or that are to 
appear in conference proceedings. A list of these reports for the period 

October 3, 1985, through March 31, 1986, is given in the Reports and Abstracts 
section which follows a brief description of the research in progress. 


Presently, ICASE is operated at NASA Langley Research Center, Hampton, VA 
under the National Aeronautics and Space Administration, NASA Contract No. 
NAS1-18107. In the past, support has been provided by NASA Contracts No. 
NAS1-17070, NAS1-17130, NAS1-15810, NAS1-16394, NAS1-14101, and NAS1-14472. 
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RESEARCH IN PROGRESS 


Saul Abarbanel 

Work of S. Abarbanel, A. Bayliss, and L. Lustman continues on the 
construction of non-reflecting boundary conditions for the compressible 
Navier-Stokes equations. The boundary conditions are obtained by linearizing 
the Navier-Stokes equations around a model mean profile. We then compute an 
approximate dispersion relation between the frequency and the waves that can 
be supported by the resulting linear system. The approximate dispersion 
relation is valid in the limit of long wavelengths. This analysis has been 
completed for flows over a flat geometry. We have verified that the new 
boundary condition can significantly accelerate the convergence to a steady 
state. A report describing this work was completed. We have completed the 
analysis for flows in a channel and are presently testing the resulting 
boundary conditions numerically. 


H. T. Banks 

The standard proofs of convergence in identification techniques for 
distributed parameter systems involves a "uniform compactness" type assumption 
on the parameter sets and their family of approximates. This compactness also 
allows one to prove "method-stability" results (e.g., continuity of parameter 
estimates with respect to observations) and has an analogue in the Tikhonov 
regularization formulation for these problems. We are currently pursuing 
investigations into the computational necessity of imposing this compactness 
through some type of functionally-constrained optimization algorithms. A 
numerical study of the algorithms including a comparison with regularization 
techniques is underway. 


H. T. Banks and K. Ito 

We have investigated the use of a hybrid method for the computation of 
feedback gains in a general class of linear quadratic Gaussian (LQG) 
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problems. The methods are easily applied to large but finite dimensional 
interconnected systems as well as approximate systems for infinite dimensional 
distributed parameter control problems. The methods entail using a 
Chandrasekhar algorithm to obtain start ups for a Newton-Kleinman algorithm 
with a variable step Smith method to solve the resulting Lyapunov equations in 
the iterative method. Numerical experiments suggest theat the method is 
superior to the standard eigenvector type techniques (Potter, Laub-Schur) in 
problems where the number of states is large compared to the number of control 
inputs . 


H. T. Banks and Robert K. Powers 

Large space structures currently proposed often contain long flexible 
members. The size of such structures makes it intractable to test full scale 
models in slewing and control maneuvers, making it necessary to have reliable 
mathematical models for computer simulation. Current research is underway in 
conjunction with J. Juang (Structures and Dynamics Division, LaRC) to model 
damping due to air during slewing maneuvers of flexible beams. Finite element 
and parameter identification techniques are being applied to models of 
flexible beams to determine the coefficients of nonlinear damping terms. 


H. T. Banks and Georg Propst 

Together with R. Silcox (Acoustics Division, LaRC) a noise suppression 
problem has been formulated as a control problem governed by the wave 
equation, and the numerical performance of approximations to this model is 
being evaluated. Several foundations of a one-dimensional version of the 
problem are being studied for computational feasibility, in addition, 
extension of the ideas to multi-dimensional domain problems are being 
considered. 
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H. T. Banks and I. G. Rosen 

We have continued our efforts on the difficult problems associated with 
the estimation of spatially varying material parameters such as stiffness and 
damping in cantilevered beams with tip mass or tip body. Results of our 
comparisons of several types of finite element state approximation schemes for 
various state formulations (e.g., weak variational, conservative) are now 
being prepared for publication. The enhancement of methods and the effects on 
"identif lability" via applied shear forces and moments at the tip have been 
investigated. 

In a related study we have developed and used similar ideas for a model of 
the "RPL experiment" (a beam with tip thrusters and a flexible hose). We have 
successfully identified physically meaningful parameters in a proof-mass 
(i.e., spring-mass-dashpot) - beam hybrid model using acceleration data from a 
series of experiments. 


Alvin Bayliss 

Work of A. Bayliss, L. Maestrello (Transonic Aerodynamics Division, LaRC) , 
P. Parikh (Vigyan Research Associates, Inc.), and E. Turkel continues on the 
numerical simulation of the active control of compressible flows by localized 
periodic surface heating. We have completed our study of the control of two- 
dimensional flows over a curved surface. The results show that the pressure 
gradient induced by the surface curvature significantly enhances the 
receptivity of the flow to surface heating. Numerical simulations demonstrate 
that unstable disturbances can be stabilized by a proper choice of the phase 
of the surface heating. A report describing this work is being prepared. We 
have begun simulating three-dimensional disturbances in a two-dimensional mean 
flow. 


Shahid H. Bokhari 

The classical (circa 200 BC) Sieve of Eratosthenes for finding prime 
numbers is one of the oldest parallel algorithms known to man. In recent 
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years it has seen much use as a benchmark algorithm for serial computers while 
its intrinsically parallel nature has gone largely unnoticed. 

We have implemented a parallel version of the Sieve on the Flex/32 
multiprocessor and are evaluating its performance. The Sieve is a good test 
of the capabilities of a multiprocessor since it requires rapid access to 
shared memory by a number of processes that are created and destroyed on the 
fly. It is also a challenge for the algorithm designer since, by the very 
nature of the algorithm, it is impossible to get good speedups without some 
form of dynamic load-balancing. We expect that this parallel version of the 
Sieve will be useful in comparing the performance of shared memory 
multiprocessors . 


Shahid H. Bokhari and Marsha Berger 

Research on a binary partitioning technique for non-uniform problems has 
been completed and a technical report prepared. We have compared the 
performance of nearest neighbor arrays and hypercubes when using our 
partitioning strategy. Some of our results are unexpected — they indicate that 
hypercubes perform better than meshes only if the problem domain has moderate 
non-uniformity. For the case of mildly distorted or heavily distorted 
domains, the performance of hypercubes is essentially equal to that of meshes. 


Shahid H. Bokhari and M. Y. Hussaini 

In collaboration with Gordon Erlebacher (High-Speed Aerodynamics Division, 
LaRC), we have embarked on a major effort to transport an existing Navier- 
Stokes code from the Cray-2 to the Flex/32. The objective is to parallelize 
the code and demonstrate speedup on the Flex/32 in anticipation of the 
availability of parallel processing on the Cray-2 and its successors. By 
carefully developing the synchronization and spawning routines needed for 
parallelism, we have succeeded in moving the code to the Flex/32 with very 
little modification. In fact, other than calls to the synchronization 
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routines and changes in the do loop ranges, there is almost no change in the 
main body of the code. This work is thus likely to have a major impact on the 
utilization of parallel machines for large fluid mechanics problems over the 
next five years. 


Shahid H. Bokhari and M. Ashraf Iqbal 

We are studying the problem of optimally assigning the modules of a 
parallel/pipelined program over the processors of a multiple computer system 
under certain restrictions on the interconnection structure of the program as 
well as the multiple computer system. We have established that for a variety 
of such programs it is possible to determine in linear time if a partition of 
the program exists in which the load on any processor is within a certain 
bound. This method, when combined with a binary search over a finite range, 
provides an approximate solution to the partitioning problem. 

Another research effort addresses the problem of finding a minimal super 
binary tree which contains a given set of given binary trees. Such a problem 
arises in the design of the ALU pipeline of the Navier-Stokes computer under 
development at Princeton University. The ALU pipeline is composed of 
individual units connected in the form of a binary tree. As each different 
application may require a distinct binary tree interconnection of the ALU 
chips, it is useful to hard wire the ALU pipeline in the form of a minimal 
super binary tree which contains all required trees. We have constructed the 
minimal supertree for a specific set of trees and are now developing an 
efficient algorithm for the general problem. 


Shahid H. Bokhari, M. Ashraf Iqbal, and Joel H. Saltz 

We have studied the problem of uniformly distributing the load of a 
parallel program over a multiprocessor system. We have compared the 
performance of four strategies for load balancing. These are (1) the optimal 
static assignment which is guaranteed to yield the best static solution, (2) 
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the static binary dissection method which is very fast but sub-optimal, (3) 
the greedy algorithm, a static fully polynomial time approximation scheme, 
which estimates the optimal solution to arbitrary accuracy, and (4) the 
predictive dynamic load balancing heuristic which uses information on the 
precedence relationships within the program and outperforms any of the static 
methods. A technical report has been prepared. 


Dennis W. Brewer 

Research is continuing on parameter estimation problems associated with 
linear evolution equations in infinite-dimensional spaces. A general 
algorithm based on quasi-linearization has been established along with its 
local convergence properties. The algorithm has been numerically tested on 
linear delay-differential equations. Numerical experiments indicate that the 
method converges rapidly when used to identify two unknown delays together 
with two or three unknown coefficients using simulated data. Future research 
will involve continuing numerical experimentation and improvements in the 
theory to accommodate a wider class of problems. 


Dennis H. Brewer and J. Steven Gibson 

We are continuing our efforts to develop robust software routines for 
identifying parameters in robotic manipulators. The programs use numerical 
integration of nonlinear differential equation models and nonlinear 
optimization algorithms. This work is in conjunction with J. Pennington and 
F. Harrison (Information Systems Division, LaRC) who have provided 
experimental data. Our previous research indicates a need to model the 
integrated electro-mechanical system to reduce parameter sensitivity. We have 
obtained data and plan to test models which include coupled physical and 
actuator dynamics. In addition, input optimization should be a fruitful area 
for future research. 
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Jan Cuny, Piyush Mehrotra, and John Van Rosendale 

BLAZE is a scientific programming language developed at ICASE over the 
past two years. BLAZE programs are mapped onto parallel architectures by a 
program transformation system. This system translates BLAZE source programs 
through a sequence of intermediate forms to Extended-BLAZE. E-BLAZE is an 
explicit-tasking language consisting of BLAZE together with a number of lower- 
level constructs for expressing parallelism. It contains constructs for 
allocating and de-allocating processes, specifying array storage patterns, 
scheduling loops, performing interprocess communication and synchronization, 
and so forth. E-BLAZE is, in effect, a high level virtual architecture for 
BLAZE. 

Though E-BLAZE is intended primarily as a compiler intermediate form, it 
may also prove useful for programmers wishing to fine-tune their algorithms 
for a particular architecture. This is especially true for non-shared memory 
architectures, such as the NCUBE or FPS hypercubes, or the Navier Stokes 
Computer being developed at Princeton University, where well developed 
software environments are currently unavailable. We are studying both the 
problem of mapping E-BLAZE to such architectures and the issue of its utility 
as an application programming language. 


Pravir Dutt and Shlomo Ta'asan 

Results have been obtained for parameter estimation problems involving 
elliptic partial differential equations using the equation error approach. 
Work is continuing on extending the method to nonelliptic problems such as the 
equations for oil recovery and the Euler-Berroulli beam equation. 


Stefan Feyock 

The previous report described the development of a Prolog escape 
predicate, which allowed convenient communication between Prolog and external 
software systems, in particular the FORTRAN-based database management system 
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RIM. The resulting integration of Prolog with a large-scale commercial 
database was used to add expert-system capabilities to the database system. A 
RIM-based expert system that diagnosed jet engine malfunctions was developed 
and run successfully. Current research is focusing on the observation that 
addition of an Artificial Intelligence (AI) capability to pre-existing 
software systems is a common requirement and that an AI capability that can be 
called in slave mode from existing software would be a valuable programming 
tool. Means of providing such interfaces between AI tools and conventional 
software are being studied. 


Dennis Gannon and Piyush Mehrotra 

The principal open questions with BLAZE revolve around the family of 
program transformations needed to map BLAZE programs to multiprocessor 
architectures. We plan to complete preliminary implementations of BLAZE on 
several types of multiprocessors during the next six months. Our first target 
architectures are the BBN Butterfly, Flex/32, and Alliant multiprocessors. We 
are also planning to implement versions of BLAZE on the IBM RP3 prototype and 
the University of Illinois Cedar computer as soon as they become available. 

Coupled with this effort to achieve preliminary implementations is a 
longer term effort to study the subtle issues involved in parallel run-time 
environments. While it is relatively easy to create naive implementations of 
BLAZE on shared memory multiprocessors, constructing implementations which 
fully exploit these architectures is difficult. For example, on the RP3 one 
could potentially exploit memory locality, caching, the fetch-add network, and 
could also allocate sub-groups of processors to cooperate on parallel tasks. 
With such architectural flexibility, deciding the best way to implement any 
given language construct can be quite complex. 

To study such issues, a graphics based program transformation system has 
been written and is currently running on Apollos and Macintoshes. This system 
will be ported to the ICASE Suns over the summer and should prove useful in 
developing program transformations suited to the NASA Flex/32 multiprocessor. 
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To some extent, this work parallels work done by D. Kuck, K. Kennedy, 
F. Allen, and others on mapping Fortran programs to parallel architectures. 
The difference is that BLAZE provides a much better starting point, since the 
semantics of BLAZE were chosen with exactly this problem in mind. A Prolog 
based expert system has been developed to select the sequence of 
transformations to be performed, and preliminary results on the efficacy of 
this approach should be available soon. 


James Geer 

During the past six months, I have been working on the following projects 
in connection with NASA Langley: 

Comparison of numerical and perturbation techniques to determine the 
aerodynamic forces on slender bodies . This study has involved the comparison 
of pressure coefficients for some two- and three-dimensional slender bodies of 
simple geometrical shape using standard numerical (panel) methods and uniform 
asymptotic methods. Some preliminary results have been reported in a recent 
ICASE report. Work is now continuing to try to determine if some of these 
results can be "pieced together" to describe the forces acting on more 
complicated geometrical shapes, such as a typical subsonic aircraft. This 
work is being done with Dr. Eddie Liu (Low Speed Aerodynamics Division, LaRC) 
and Prof. Lu Ting of New York University. 

Reduced basis method extensions and applications . This work involves a 
study of the reduced basis method (a semi-numerical, semi-analytical method 
which has been used to solve a variety of nonlinear boundary value problems in 
elasticity) to determine how the method might be applied to some fluid dynamic 
and scattering problems. Several model problems are being studied both to 
understand the method better and to see what type of problems the method might 
be applied to. At present, these model problems include some two point 
boundary value problems with boundary layers and some simple exterior 
scattering problems. This work is being done with Dr. Carl Andersen of the 
College of William and Mary. 
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Nonlinear oscillations using symbolic computation . This study involves 
the use of the symbolic computation system MACSYMA and perturbation methods to 
investigate some free and forced nonlinear oscillations. The method of 
multiple time scales has been implemented using MACSYMA and applied to the 
problem of determining the transient responses of the van der Pol and Duffing 
oscillators when the applied force consists of a sum of periodic terms with 
different frequencies. The conditions under which the system will experience 
frequency entrainment are being investigated. This work is being done with 
Dr. Carl Andersen of the College of William and Mary. 

Integral equations of the first kind . Investigations are underway 
concerning some fundamental properties (both from theoretical and 
computational points of view) of a class of "almost" singular integral 
equations of the first kind which are useful in representing the solution to 
certain elliptic exterior boundary value problems. These equations typically 
have the property that the domain of integration R is a proper subset of the 
domain of validity D of the equation. Special consideration is being given 
to the idea of analytically continuing the solution into the domain D. 
Applications of the results will be useful in several two- and three- 
dimensional problems involving slender or thin bodies. 


J. S. Gibson and I. G. Rosen 

We are developing numerical approximation methods for the computation of 
optimal feedback control and observer gains for infinite dimensional, 
discrete-time, linear-quadratic regulator problems. In particular, we are 
considering systems in which the state transition operators are given in terms 
of linear semigroups of operators on infinite dimensional Hilbert spaces and 
the input and output operators are possibly unbounded. Our investigation is 
to include both theoretical (convergence arguments, etc.) and numerical 
components. We intend to apply our theory to a variety of test examples 
including systems describing transport phenomena (e.g., diffusion, heat 
conduction, etc.) and the vibration of flexible structures. 
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We are also studying numerical approximation methods for the estimation or 
identification of spatially varying coefficients in wave and beam equations 
based upon model data (e.g., natural frequencies and mode shapes). 
Specifically, we are looking at schemes which use spline-based finite element 
methods to discretize the infinite dimensional inverse problem. A scheme 
which uses modal data rather than time histories of displacement, velocity, 
strain, etc. to estimate unknown parameters has the potential to be much 
faster since the repeated integration of the initial-boundary value problem 
which governs the underlying dynamics (or more precisely, a high order finite 
dimensional approximation to it) is not required when the optimization is 
carried out. 


David Gottlieb 

Research continued on the problem of the behavior of high order schemes 
applied to discontinuities. It is clear that the notion of formal accuracy is 
not valid and, in fact, loss of accuracy is observed. During the last few 
years, we have realized that high order schemes do contain more information 
about the solution if they can be made stable. We are still working on the 
general theory, and we are making progress in applying the notion of 
"information content" to high order schemes. Joint work with D. Fishelov 
(University of California, Berkeley) demonstrates that accurate results can be 
recovered from spectral simulations of the small disturbance equation. 

Together with M. Salas (Theoretical Aerodynamics Division, LaRC) research 
is in progress concerning the vortex breakdown problem. We study properties 
of solutions of the ODE model developed by Salas. 

Together with T. Zang and G. Erlebacher (Computational Methods Division, 
LaRC) we study splitting methods for time dependent flows with small 
variations in the sound speed. The idea is to remove the stiffness of the 
system by splitting out the average (constant) sound speed. 
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Chester G. Grosch 

A three-dimensional, time dependent Navier-Stokes code in vorticity and 
velocity variables has been developed in collaboration with T. Gatski (Viscous 
Flow Division, LaRC) and M. E. Rose. This code is being modified in order to 
increase its speed and efficiency and is being tested against several exact 
analytic solutions. The initial calculations for the flat plate receptivity 
problem have been completed. The results have been analyzed, written up, and 
submitted for publication. 

The study, in collaboration with R. Fatoohi, of parallel algorithms for 
elliptic problems is continuing. The four color relaxation scheme is now 
working on the MPP at NASA Goddard. This has been quite successful for 
problems which fit on the array; the average processing rate for these 
problems is 175 MFLOPS. A modification of the algorithm for problems which 
are much larger than the array has been developed. This will be tested in the 
near future. Similarly, a modified algorithm has been developed for the MIMD 
FLEX/32 and is being tested. 


Subramaniya 1. Hariharan 

Currently, we are making numerical simulation of atmospheric acoustic 
problems that were theoretically investigated during the last reporting 
period. The computer code for this problem is stable and consistent with the 
theoretical boundary conditions that we proposed. Results will be tested 
against an experiment that is planned in the fall of 1986 at Langley. 

The major project that we are undertaking here is to study the effect of 
acoustics on flows past blunt bodies. Even the steady compressible viscous 
flow problem in the situation of "leading edge" is quite complicated. 
Acoustics is an unsteady phenomenon. Thus, as a first attempt we are 
investigating efficient methods, such as optimal time splitting, to obtain 
simulations of flow past a flat plate including the leading edge. Then the 
acoustic disturbances will be introduced in the flow. In this effort, well- 
posedness of the problem and appropriate outflow boundary conditions with and 
without acoustics will be examined. 
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Yousuf f Hussain! 

In a continuing program of research in hydrodynamic stability and 
transition, refined simulations of certain basic problems were carried out 
with a view to unravelling the underlying instability mechanisms. In 
particular, the simulation of the turbulent Rayleigh-Benard problem on a 124 x 
64 x 64 grid yielded results in good agreement with experiments; these are 
reported in ICASE Report No. 86-6. Statistical analysis of the simulation 
data is continuing to assess the effect of buoyancy on subgrid scale 
modelling, etc. Another fundamental stability problem under study is the 
Taylor-Couette flow in a very short annulus allowing for an exchange between 
two-cell and single-cell modes. Of immediate interest is the temporal 
evolution of the order parameter which quantifies the symmetry-breaking 
bifurcation of a symmetric two-cell flow into an asymmetric single-cell 
flow. The long term objective is the simulation of the three-dimensional 
processes leading to transition. 

Our investigation of the receptivity problem is continuing. We have 
constructed an accurate and efficient spectral method of solving the three- 
dimensional compressible Navier-Stokes equations relevant to certain physical 
experiments. The purpose is to assess the effect of Mach number on the 
instability mechanisms identified in low-speed flows. 

D. Bushnell (Viscous Flow Division, LaRC) proposes a method of suppressing 
the lift-induced drag. It consists of essentially pushing down the front 
stagnation point and pushing the rear stagnation point up on the upper 
surface. This concept is being tested in the framework of the two-dimensional 
compressible Navier-Stokes simulation. 

This program of research is being carried put in collaboration with G. 
Erlebacher, T. A. Zang (High-Speed Aerodynamics Division, LaRC), and C. L. 
Streett (Transonic Aerodynamics Division, LaRC). 


Kazufumi Ito and Robert K. Powers 

Efforts are continuing in the application of Chandrasekhar-type algorithms 
to infinite dimensional linear quadratic regulator (LQR) problems that have 
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unbounded inputs and outputs (i.e., boundary control and observation). 
Approximation of the associated gain operators via the Chandrasekhar equations 
has been shown to be more computationally efficient than approximation via the 
Riccati equation for many problems with bounded inputs and outputs. The 
features of the Chandrasekhar equations that make them computationally 
efficient are present in many problems involving boundary controls and 
observations. Numerical techniques for solving such problems are being 
sought. Research on application of these ideas to LQR problems with unbounded 
input and output operators is also being pursued. 


Tom L. Jackson 

An earlier paper titled "Shock- Induced Thermal Runaway" (SIAM J. Appl. 
Math., 1985) describes the initial stage in the birth of a detonation wave. 
The ignition stimulus was modelled by impulsively moving a piston into a 
combustible gas confined within a tube at constant speed. This creates a 
shock wave which runs ahead of the piston and ignites the gas at the piston 
face. The creation of the hot spot and its growth into a full-fledged 
explosion in the vicinity of the piston face was described. We are currently 
investigating techniques to develop the solution further in time, describing 
the complete history of transition to ignition. This work is being done in 
collaboration with A. K. Kapila. 


Tom L. Jackson and Pravir Dutt 

A mathematical model is being constructed to examine the combustion 
process that occurs when two gases, initially separated by a finite length 
splitter plate, come into contact and mix at the trailing edge. For large 
Reynolds numbers, the trailing edge region can be represented analytically by 
the so-called "triple deck" model of Stewartson and Messiter. It is found 
that mixing and ignition takes place only in the lower deck and can be 
described by the standard boundary layer equations subject to unusual boundary 
conditions. A code is currently being constructed to solve this set of 
equations. This work is being done in collaboration with M. Y. Hussaini. 
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Harry Jordan 

The parallel programming environment known as the Force and described in 
ICASE Report 85-45 has been implemented under the Concurrent Fortran system on 
the FLEX/32 computer. Programs previously run on the HEP multiprocessor have 
been used to test the implementation of the system. Work is currently 
underway to parallelize some of the Fortran routines in the NICE/SPAR system 
using the parallel extensions embodied in the Force. Of particular interest 
is whether extensions to the Force system are indicated as a result either of 
the FLEX/32 hardware and performance parameters or of the type of algorithms 
appearing in NICE/SPAR. 


William D. Lakin 

Instability of a two-dimensional boundary layer flow through the 
interaction of inviscid- and viscous-type disturbances . The basic flow in 
this investigation is a two-dimensional, parallel boundary layer profile with 
non-zero mean velocity components in the streamwise and spanwise directions. 
When a small amplitude, three-dimensional, inviscid-type disturbance is 
imposed on this mean flow, the equations governing the amplitudes of the 
velocity and pressure perturbations can be reduced to a Rayleigh equation 
involving two wavenumbers. The present work considers possible instabilities 
due to the interaction of these inviscid-type disturbances with viscous-type 
disturbances. 

The initial focus of this work has been on understanding the nature of the 
inviscid disturbances. Frobenius-type solutions of the Rayleigh equation have 
been derived. The asymptotic behavior of the inviscid disturbances has also 
been obtained making possible determination of the appropriate combination of 
inviscid perturbations which remains bounded in the free stream. The 
possibility of resonant inviscid triads is currently being investigated. This 
work is being done in collaboration with P. Hall and M. Y. Hussaini. 

Determination of realistic modal damping coefficients for a rotor-airframe 
helicopter model . Theoretical investigations of the forced response of a 
helicopter, modeled as a rotor-airframe structure, typically are based on a 
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finite-element model in which the modal damping coefficients are taken to be 
uniform for all modes. Recent comparison of theoretical predictions with 
experimental data shows that uniform modal damping produces poor correlation 
between theory and experiment. The present work investigates the feasibility 
of determining realistic, non-uniform damping coefficients which will bring 
about close correlation with experimental response curves in frequency bands 
away from resonant peaks. This work is being done jointly with Ray Kvaternik 
(Loads and Aeroelasticity Division, LaRC). 


Liviu Lust man 

The steady state elasticity problem for composite materials has been 
solved for the simple brick geometry. The compact scheme equations are solved 
by conjugate gradient. The numerical results compare favorably with 
calculations performed on the same problem by other algorithms. 

Work is in progress on the more realistic problem of a plate with 
stiffeners, which can be reduced to the previous one — the domain considered 
may be decomposed into two bricks with a common interface. The boundary data 
for this problem are to be supplied from a two-dimensional finite element 
code, and the results — a three-dimensional analysis — may validate the two- 
dimensional model and add details to it. This work has been performed in 
collaboration with M. E. Rose and S. Ta'asan. 

Joint research with D. Gottlieb and E. Tadmor has continued on spectral 
collocation methods applied to initial-boundary value problems for hyperbolic 
systems. We can prove convergence of the spectral approximations to the 
solutions of the differential equations. 

The study of optimal numerical boundary conditions for steady state 
Navier-Stokes is continuing, in collaboration with S. Abarbanel and A. 
Bayliss. We have identified the parameters for channel flow, using a model 
O.D.E. eigenvalue problem. Test calculations are in progress for the full 
Navier-Stokes code with these novel boundary conditions. 
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Nessan Mac Giolla Mhulris 

Abrupt and drastic changes of structure can often occur in vortex flows. 
Theoretical investigations of such phenonema have focused on axisymmetric 
processes. Nevertheless, it has been experimentally verified that all 
transitions occurring in concentrated vortex flows are three-dimensional in 
nature. We will try and simulate numerically the transition sequence that has 
been observed experimentally in the swirling flow of an incompressible fluid 
in a pipe. Two fundamental parameters govern such a flow, namely, the amount 
of swirl imparted to the flow and its Reynolds number. As these parameters 
are altered from low values, flow visualization studies indicate that the 
axisymmetry of the initial laminar flow is broken and, while the ultimate form 
of transition to the "bubble" type of "vortex breakdown" (encountered when the 
governing parameters reach their largest values) appears to the eye to be 
axisymmetric, it is not, nor are any of the other structures that appear at 
lower values of the parameters. 

Our investigation will focus on the three-dimensional aspects of these 
flows. The numerical method employed in the study will be a spectral one that 
is ideally suited for the simulation of incompressible flows. The velocity of 
the fluid is expanded in a series of divergence free vectors which ensures 
that the continuity equation for an incompressible fluid is satisfied exactly 
and the pressure term is eliminated from the Navier-Stokes equation by 
projecting it onto the space of inviscid vector fields. By choosing 
appropriate polynomials for the components of the divergence free vectors, we 
can achieve "exponential accuracy" for the components of the velocity. Such 
expansions have been used successfully in the past to solve certain stability 
problems and also for the simulation of fluid flows such as the well known 
Taylor-Couette flow. 


Piyush Mehrotra and John Van Rosendale 

BLAZE is a scientific programming language designed at ICASE during the 
last two years. It has a Pascal-style syntax, but employs functional 
procedure calls, making it relatively easy to restructure BLAZE programs for 
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execution on multiprocessor architectures (The BLAZE Language, ICASE Report 
No. 85-29). There are several advantages to this approach. It frees 
programmers from the necessity of dealing with the complex parallel run-time 
environment implicit in the use of explicit-tasking languages, such as Ada or 
Occam, and also allows the expression of parallel scientific programs at a 
higher and more natural level. 

A BLAZE compiler which produces sequential code for Vaxes is currently 
available, and a number of parallel implementations of BLAZE are being 
developed. We are currently in the process of defining Version 2 of BLAZE. 
This will differ from the current language primarily through the introduction 
of abstractions and compilation modules. Using abstractions, one will be able 
to construct large objects, such as files, which have their state retained 
between procedure calls making it easier to write large scientific programs in 
BLAZE. 


Vi jay K. Naik 

In the past two years, a variety of general purposes MIMD architectures 
has appeared as commercial products. Many algorithms suitable for parallel 
processing have been proposed. The development, on parallel machines, of 
application codes for solving problems of practical importance are also 
reported. In spite of these, several difficulties exist in implementing 
nontrivial problems efficiently on the multiprocessor systems. These 
difficulties are primarily because of the lack of methodologies and tools that 
can automate the process of matching the rates and the amounts of data that 
must be transferred among the various computational tasks with the structure 
and the capacity of the underlying interconnection network. 

The study undertaken here is aimed at developing a methodology that will 
allow one to characterize MIMD architectures and will quantify the 
architecture and algorithm dependent parameters that contribute to the 
communication delays. Various existing, as well as proposed, architectures 
and important algorithms are currently under investigation. The outcome of 
this research will help towards the automation process of efficient parallel 
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implementations so that the user need not know the detailed anatomy of the 
underlying architecture. The research currently in progress with Shlomo 
Ta'asan (described below) addresses some of the issues involved in solving 
scientific problems on various currently existing architectures. 


Vi jay K. Naik and Shlomo Ta'asan 

Algorithms based on multigrid techniques are known to be among the most 
efficient methods for solving partial differential equations. On sequential 
machines, a wide class of problems discretized on a grid with n points 
requires only 0(n) arithmetic operations to obtain a solution to within the 
truncation error of the discretization. This makes the multigrid methods 
optimal on the sequential machines. This optimal property and the fact that 
many aspects of the multigrid algorithms are highly parallelizable , make these 
methods attractive for implementation on multi-processor systems. Several 
researchers have studied the effects of parallelizing the intra-grid and 
inter-grid operations. Some researchers have proposed special purpose 
architectures that allow one to map some important problems directly onto the 
hardware and thus reduce the communication delay effects. But none of these 
studies have conclusively shown that the optimality observed on the sequential 
machines is conserved in time and the number of processors needed when these 
methods are implemented on multiprocessors, nor is there any clear 
understanding about the communication and the synchronization costs involved. 

The work undertaken here is to explore the above mentioned questions 
further for shared and non-shared memory architectures. The effort is aimed 
at studying the effects of various communication strategies and also at 
understanding the effects of various implementations on the convergence 
rates. We are also studying the advantages and disadvantages of using a large 
number of simple processors and those of using a small number of powerful 
processors. The model problem considered here is that of solving 2-D 
incompressible Navier-Stokes equations in vorticity-stream function 
formulation. Experiments are being carried out on commercially available MIMD 
machines such as the Intel Hyper-cube, Sequent, Encore, and Flex/32. In 
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addition to these studies, investigation is underway towards the development 
of new algorithms that possess some of the numerical properties of multigrid 
methods but do not require tight synchronizations. 


David M. Nicol 

A large class of parallel computations have run-time behavior which may 
vary in time. Examples of such computations include many numerical methods 
for solution of scientific problems and parallel simulations. Because of this 
stochastic variance, a good partitioning of the workload between processors 
can drift out of balance, leading to degraded performance. This performance 
degradation gives rise to two fundamental problems: (1) how to model this 
phenomenon, and (2) how to detect and respond to performance degradation due 
to stochastic variation. Future efforts will be focused on both of these 
issues. 


Merrell L. Patrick 

Computing systems of the future for large scale scientific and engineering 
computations will consist of many processors and memories working in parallel 
to solve a single problem. Efficient use and ease of use of these future 
systems are two basic problems facing potential users. The development of new 
parallel algorithms, parallel execution time models which predict the 
performance of the algorithms on different parallel architectures, and 
parallel programming environments for expressing the algorithms are essential 
for progress toward the solution of these problems. 

We are developing parallel algorithms and corresponding parallel execution 
time models for solving the generalized eigenvalue problem. Results from the 
models will be used to predict the performance of the algorithms on message 
passing and shared memory MIMD machines. The parallel algorithms will be 
integrated into a structural analysis system, NICE/SPAR, running on the 
FLEX/32 as a program module used in the dynamic analysis of a structure in 
motion. 
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The algorithms developed will be used as a basis to continue work with 
T. Pratt and H. Jordan on defining the essential ingredients of parallel 
programming languages for scientific computations. Different programming 
paradigms will be used to program these algorithms on the FLEX/32 system. 

Recent work with L. Adams and D. Reed determined optimal stencil/partition 
pairs for solving elliptic partial differential equations on MIMD machines. 
The optimality condition was to maximize the computation to communication 
ratio. Parallel execution time models were developed for message passing and 
shared memory architectures to predict the impact of the architecture on the 
choice of stencil/partition pairs. The results of the models will be 
validated using various available machines such as the FLEX/32, Sequent 
Balance 8000, and the Connection Machine. 

A new study was begun with Saltz and Voigt which attempts to answer the 
question as to whether SIMD architectures are sufficient for carrying out most 
computations arising in the solution of partial differential equations. A 
thorough review of fundamental algorithms used in solving partial differential 
equations will be conducted in an attempt to identify those which require MIMD 
architectures . 


Doug Peterson 

The VAX and Ridge computer systems, as well as their installation 
configurations, both hardware and software have been learned. The Talaris 
laser printer has been brought on-line, and has been used with ditroff, pic, 
and System V graphics. The VAX has provided user support to the NAS project 
by sharing the tape drive with the Iris workstation as a network resource and 
as an interactive gateway to the computing resources at Ames Research Center 
through the micom switching system. 

To better allocate ICASE computing resources, a patch panel was designed 
and installed in the VAX room to provide quick and reliable interconnections 
between user terminals, the VAX, and micom. The system has been documented as 
both a set of floor plan drawings and as a table of connections. To provide 
for future updates, the drawings are done with the Macintosh using MacDraw 
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software, and the tables are kept on the VAX. A similar set of floor plans 
has been developed for the ICASE office staff to aid in developing seating 
plans. 

I have continued to provide consultation and support for other UNIX system 
administrators, both at Langley and other locations. User support for Langley 
is being absorbed by the Analysis and Computation Division, and I expect to 
reduce that activity at ICASE. However, the locally prepared manual 
"Introduction to Computing at ICASE" is being updated. I have offered to 
participate in the addition of UNIX systems to the LarcNet. 

Plans have been made to replace the VAX with a system of Sun 3 diskless 
workstations and a file server, connected via ethernet. Delivery is expected 
by May 9, 1986, with installation on or about May 15, 1986. 


Terrence W. Pratt 

The PISCES parallel programming environment allows the scientific 
programmer to write parallel programs that are intended to run on various 
parallel architectures. Programming is in PISCES Fortran (an extension of 
Fortran 77) under the Unix operating system. Several "granularities" of 
parallel operations are provided. 

The PISCES 1 environment has been implemented on the ICASE VAX 750 and on 
a network of Apollo Domain workstations at the University of Virginia. In 
PISCES 1, a program is organized in clusters of tasks that communicate by 
passing messages. 

The PISCES 2 environment is being designed and implemented for the FLEX/32 
at NASA Langley. PISCES 2 extends PISCES 1 to include the 'force' constructs 
developed by H. Jordan which allow the programmer to make effective use of the 
shared memory of the FLEX/32. PISCES 2 also allows the programmer to control 
the mapping from PISCES clusters to the hardware processors. 
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Georg Propst 

Approximation schemes for hereditary differential equations are studied. 
In order to prove applicability of the approximations to control problems, 
strong convergence of the adjoint systems is required. One subject of the 
current work is the investigation of the adjoints of a class of spline 
approximations. Another subject is a scheme that uses discontinuous piecewise 
linear functions and has favorable stability properties. If the cost 
functional ranges over an infinite time horizon, the proof of convergence of 
the approximating controls relies on the nonsingularity of the weight on the 
state trajectory. The necessity of this assumption is investigated. 


Joel H. Saltz, Vijay Naik, and David M. Nicol 

The efficient implementation of algorithms on multiprocessor machines 
requires that the effects of communication delays be minimized. The effects 
of these delays on the performance of a model problem on a hypercube 
multiprocessor architecture is investigated, and methods are developed for 
increasing algorithm efficiency. The model problem under investigation is the 
solution by red-black Successive Over Relaxation of the heat equation; most of 
the techniques to be described here also apply equally well to the solution of 
elliptic partial differential equations by red-black or multicolor SOR 
methods. 

A complete report identifies methods for reducing communication traffic 
and overhead on a multiprocessor and presents the results of testing these 
methods on the Intel iPSC Hypercube. We examine methods for partitioning a 
problem's domain across processors, for reducing communication traffic during 
a global convergence check, for reducing the number of global convergence 
checks employed during an iteration, and for concurrently iterating on 
multiple time-steps in a time dependent problem. Our empirical results show 
that use of these methods can markedly reduce a numerical problem's execution 
time. 
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Joel H. Saltz and David M. Nicol 

We consider the tradeoffs between the benefits and costs of balancing load 
in multiprocessor solutions of computational problems that describe the 
evolution of a physical system over the course of time. The type of 
architecture being assumed is one in which each processor has a local memory 
and that the processors pass messages through either a communications network 
or a shared memory. The responsibility for computations pertaining to a given 
portion of problem domain is assigned to each processor. The data 
dependencies inherent in the solution of a time dependent problem cause the 
most heavily loaded processor to limit the rate at which the problem is 
solved. Frequently, the load in the portion of the problem assigned to each 
processor will vary, and without redistribution of load, system performance 
declines. While performance can be improved by a redistribution of load, this 
redistribution itself exacts a sometimes large delay cost. 

Two probabilistic load models are proposed and used to evaluate policies 
for deciding when load balancing should be performed. Through use of these 
models, it is possible to characterize the relative performance of a number of 
different policies designed to determine when load should be balanced. This 
performance characterization is carried out through the derivation of analytic 
expressions and through simulation. Estimates of how the cost of balancing 
load interact with model parameters describing the number and activity of 
processors are obtained. A description of this work entitled: "A Statistical 
Methodology for the Control of Dynamic Load Balancing" is to be published in 
the SIAM publication Parallel Processing and Medium Scale Multiprocessors. 


Joel H. Saltz and Hillel Tal-Ezer 

The use of spectral methods in space to solve time dependent problems 
makes the investigation of efficient high accuracy schemes to march the 
solution forward in time quite attractive. Richardson extrapolation of 
leapfrog and Crank Nicholson time integration methods for the solution of 
these partial differential equations can produce efficient high order schemes 
and can be used in a number of different ways to facilitate efficient 
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solutions of spectral problems on a variety of multiprocessor architectures. 
Our investigation takes into account both numerical convergence as well as 
algorithm performance on message passing and shared memory multiprocessors. 

A detailed analysis is carried out of a model hyperbolic problem solved by 
the pseudospectral Fourier method in space and extrapolated Leap-Frog scheme 
in time. Estimates of parallel efficiency are calculated from the analysis of 
several ways of mapping the model problem onto a range of message passing and 
shared memory architectures. Given a desired error tolerance for time 
integration, taking into account the multiprocessor machine architecture, and 
the number of mesh points used to solve the model problem, we have developed a 
model that predicts the optimal number of extrapolations that should be 
carried out. Experimental results on the convergence and stability of the 
extrapolated algorithms have been encouraging, and an experimental analysis of 
the performance of these algorithms on the Flex/32 multiprocessor is being 
carried out. 


Ke-Gang Shih 

Study of a third order nonlinear ordinary differential equation describing 
laminar flow through porous channels with injection is reported in ICASE 
Report 86-17. 

Problems under consideration are nonlinear singular perturbation problems 
of chemical kinetics, enzyme kinetics, etc. We hope to get good estimates to 
help the numerical treatment of such equations. 


Shlomo Ta'asan 

Multigrid treatment of highly oscillatory problems, such as exterior 
problems in acoustics, is the subject of another study. The main idea here is 
to use a ray representation of the solution (similar to geometrical optics) on 
coarse grids, which allows a natural treatment of radiation boundary 
conditions as well as other advantages. This work is in collaboration with A. 
Brandt (Weizmann Institute). 
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Another direction of multigrid development is for system identification. 
Here the problem is to estimate certain coefficients in an elliptic equation 
for which a right-hand side and a solution are given. Estimation of the 
coefficient is needed under some minimization constraint. 


Eitan Tadmor 

Work on the approximate solution of the gas dynamics and other 
conservative systems is underway. We have implemented our proposed second 
order accurate nonoscillatory central differencing code obtaining very 
reliable results. Research is continuing on the application of entropy stable 
approximations to those systems based on the entropy conservative finite 
element method devised earlier. In particular, ideas of compensated 
compactness are used in order to control oscillations. A minimum 
(therodynamic) entropy principle was derived for approximate Euler solvers. 

We have developed a stability theory for (pseudo-) spectral approximations 
for the approximate solution to hyperbolic mixed initial-boundary value 
systems with D. Gottlieb and L. Lustman. Our stability normal-mode analysis 
is demonstrated for Chebyshev and Legendre collocation methods, and 
consequently convergence follows. Together with M. Goldberg (Israel Institute 
of Technology), we have devised new, easily checkable sufficient stability 
criteria for difference approximations to such mixed systems, extending 
previous work. 


Hillel Tal-Ezer 

In previous ICASE reports (84-8, 85-9), we have described a spectral 
algorithm (both in time and space) for solving partial differential 
equations. It has been shown to be an efficient tool for problems when the 
eigenvalues of the spatial operator have a specific orientation: either close 
to the imaginary axis or close to the negative real axis. Present research is 
focused on implementing the approach discussed in 84-8 for these kind of 
problems . 
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Eli TUrkel 

Work is continuing with C. Swanson (Transonic Aerodynamics Division, LaRC) 
on the properties of residual smoothing. The effects of the Navier-Stokes 
terms were analyzed. In addition, we studied the application of a Crocco type 
scheme to accelerate the convergence of the Runge-Kutta algorithm for viscous 
flow. Computations demonstrated that increased time steps could be realized 
for both laminar and turbulent problems. An ICASE report was completed on 
this. Further work is continuing on the analysis of residual smoothing for 
nonuniform meshes. Work is also continuing with V. Vatsa (Transonic 
Aerodynamics Division, LaRC) on the three dimensional viscous code. An ICASE 
report was completed with A. Kumar (High-Speed Aerodynamics Division, LaRC) 
and A. Moitra (Systems and Applied Sciences) on the extension of Runge-Kutta 
schemes to nozzle geometries for both subsonic and supersonic flows. 

A joint project with S. Yaniv and U. Landau (Tel Aviv University) was 
begun to study the effect of nonuniform meshes on the accuracy and convergence 
rates of schemes. It was found that for some finite volume schemes that the 
method could even be inconsistent for arbitrary meshes. In general, when the 
meshes are not sufficiently smooth, there is a loss of accuracy in physical 
space even though the accuracy in computational space is not affected. Cell 
centered schemes are more sensitive to nonsmooth grids than nodal schemes. A 
way of increasing the accuracy by accounting for the nonuniform mesh while 
keeping conservation was proposed. Two ICASE reports were completed. 

The study of preconditionings for slow flow was continued. Comparisons 
with incompressible preconditionings were also studied. Several different 
preconditionings for both compressible and incompressible steady flow were 
compared. An ICASE report was completed. 

Joint work is. continuing with A. Bayliss, L. Maestrello (Transonic 
Aerodynamics Division, LaRC) and P. Parikh (Vigyan Research Associates). This 
is a project to analyze time dependent compressible laminar viscous boundary 
layers by numerically solving the complete Navier-Stokes equations. A Fourth 
order, in space, code is used to resolve the unstable modes. Extensions were 
completed to treat flows over several ramp-like geometries. Unstable modes 
were found for several of these cases when the inflow was perturbed by a time 
periodic disturbance. Both passive and active controls could be found to 
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stabilize the flow over the variable geometry domains. Several ICASE reports 
appeared describing the method and the computational results of the control. 
The three dimensional code has been debugged and preliminary results have been 
obtained. 



REPORTS AND ABSTRACTS 


Jordan, H. F.: Vaaaltct computation with the fa on.ee . ICASE Report No. 85-45 , 

October 4, 1985, 12 pages. Submitted to Proc. of 1986 International 

Conference on Parallel Processing. 

A methodology, called the force, supports the construction of programs 
to be executed in parallel by a force of processes. The number of 

processes in the force is unspecified but potentially very large. The 
force idea is embodied in a set of macros which produce multiprocessor 
Fortran code and has been studied on two shared memory multiprocessors of 
fairly different character. The method has simplified the writing of 

highly parallel programs within a limited class of parallel algorithms and 
is being extended to cover a broader class. This paper deals with the 
individual parallel constructs which comprise the force methodology. Of 
central concern are their semantics, implementation on different 
architectures, and performance implications. 


Hall, P.: I n&tabitity o fa timc-pcniodic fatom . ICASE Report No. 85-46 , 

Ocotber 10, 1985, 28 pages. Submitted to Proc. of ICASE/NASA Workshop on 
Stability of Time Dependent and Spatially Varying Flows, August 19-20, 
1985, Hampton, VA. 

The instabilities of some spatially and/or time-periodic flows are 
discussed; in particular, flows with curved streamlines which can support 
Taylor-Gortler vortices are described in detail. The simplest flow where 
this type of instability can occur is that due to the torsional 
oscillations of an infinitely long circular cylinder. For more complicated 
spatially varying time-periodic flows, a similar type of instability can 
occur and is spatially localized near the most unstable positions. When 
nonlinear effects are considered, it is found that the instability modifies 
the steady streaming boundary layer induced by the oscillatory motion. It 
is shown that a rapidly rotating cylinder in a uniform flow is susceptible 
to a related type of instability; the appropriate stability equations are 
shown to be identical to those which govern the instability of a Boussinesq 
fluid of Prandtl number unity heated time periodically from below. 


Hall, P.: Tayloa-Goftilefi initabiiiti&i, ofa Totlmicn-Schtichting wav&i and 

othe,K fatom Qoocnncd by the intcKactioc boundary Zaye/i equations . ICASE 
Report No. 85-47 , October 11, 1985, 29 pages. To appear in J. Fluid Mech. 

The Taylor-Gortler vortex instability equations are formulated for 
steady and unsteady interacting boundary layer flows of the type which 
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arise in triple-deck theory. The effective Gortler number is shown to be a 
function of the wall shape in the boundary layer, and the possibility of 
both steady and unsteady Taylor-Gortler modes exists. As an example the 
steady flow in a symmetrically constricted channel is considered, and it is 
shown that unstable Gortler vortices exist before the boundary layers at 
the wall develop the Goldstein singularity discussed by Smith and Daniels 
(1981). As an example of an unsteady spatially varying basic state, we 
also consider the instability of high frequency large amplitude Tollmien- 
Schlichting waves in a curved channel. It is shown that they are unstable 
in the first 'Stokes layer stage' of the hierarchy of nonlinear states 
discussed by Smith and Burggraf (1985). The Tollmien-Schlichting waves are 
shown to be unstable in the presence of both convex and concave curvature. 


Alaghband, G. and H. F. Jordan: MudLti.pfiOC.Ui> OH bpoJ\AZ L/U decomposition with, 

contfiottzd i-LLL-in. ICASE Report No. 85-48 , October 18, 1985, 38 pages. 
Submitted to IEEE Trans. Comput. 

During L/U decomposition of a sparse matrix, it is possible to perform 
computation on many diagonal elements simultaneously. Pivots that can be 
processed in parallel are related by a compatibility relation and are 

grouped in a compatible set. The collection of all maximal compatibles 
yields different maximum sized sets of pivots that can be processed in 
parallel. Generation of the maximal compatibles is based on the 

information obtained from an incompatible table. This table provides 
information about pairs of incompatible pivots. In this paper, generation 
of the maximal compatibles of pivot elements for a class of small sparse 
matrices is studied first. The algorithm involves a binary tree search and 
has a complexity exponential in the order of the matrix. Different 
strategies for selection of a set of compatible pivots based on the 
Markowitz criterion are investigated. The competing issues of parallelism 
and fill-in generation are studied and results are provided. A technique 
for obtaining an ordered compatible set directly from the ordered 
incompatible table is given. This technique generates a set of compatible 
pivots with the property of generating few fills. A new huerlstic 
algorithm is then proposed that combines the idea of an ordered compatible 
set with a limited binary tree search to generate several sets of 

compatible pivots in linear time. Finally, an elimination set to reduce 
the matrix is selected. Parameters are suggested to obtain a balance 

between parallelism and fill-ins. Results of applying the proposed 
algorithms on several large application matrices are presented and 
analyzed. 


Speziale, C. G.: Szcond-ofidc.fi cLosafiz models ion rotating tufibulent filows . 

ICASE Report No. 85-49 , October 21, 1985, 24 pages. Submitted to Phys. 
Fluids. 
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The physical properties of the commonly used second-order closure 
models are examined theoretically for rotating turbulent flows. 
Comparisons are made with results which are a rigorous consequence of the 
Navier-Stokes equations for the problem of fully-developed turbulent 
channel flow in a rapidly rotating framework. It is demonstrated that all 
existing second-order closures yield spurious physical results for this 
test problem of rotating channel flow. In fact, the results obtained are 
shown to be substantially more unphysical than those obtained from the 
simpler K-C and K-& models. Modifications in the basic structure of these 
second-order closure models are proposed which can alleviate this problem. 


Kreiss, g. and H.-o. Kreiss: Convergence to steady i>tate o( ) solutions oi 
Borgeri 1 liquation. ICASE Report No. 85-50 , December 10, 1985, 30 pages. 

Consider the initial-boundary value problem for Burgers' equation. It 
is shown that its solution converge for t, to a unique steady state. The 
speed of the convergence depends on the boundary conditions and can be 
exponentially slow. Methods to speed up the rate of convergence are also 
discussed. 


Tadmor, E.: The. numerical vti>coi>iAy oi entropy stable 6cheme& ion iy&temi oi 

conservation tarn. I. ICASE Report No. 85-51 , November 8, 1985, 30 

pages. Submitted to Math. Comp. 

Discrete approximations to hyperbolic systems of conservation laws are 
studied. We quantify the amount of numerical viscosity present in such 
schemes, and relate it to their entropy stability by means of comparison . 
To this end, conservative schemes which are also entropy conservative are 
constructed. These entropy conservative schemes enjoy second-order 
accuracy; moreover, they admit a particular interpretation within the 
finite-element framework, and hence can be formulated on various mesh 
configurations. We then show that conservative schemes are entropy stable, 
if and — for three-point schemes — only if, they contain more viscosity than 
that present in the mentioned above entropy conservative ones. 


Maestrello, L. , P. Parikh, A. Bayliss, and E. Turkel: kctive control oi 

compressible ilom on a curved suriace. ICASE Report No. 85-52 , 
November 14, 1985, 28 pages. 

We consider the effect of localized, time-periodic surface heating and 
cooling over a curved surface. This is a mechanism for the active control 
of unstable disturbances by phase cancellation and reinforcement. It is 
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shown that the pressure gradient induced by the curvature significantly 
enhances the effectiveness of this form of active control. In particular, 
by appropriate choice of phase, active surface heating can completely 
stabilize an unstable wave. 


Harabetian, E.: Matched asymptotic, expansions to similarity solutions of 

shock diffraction. ICASE Report No. 85-53 , November 18, 1985, 32 pages. 
Submitted to Comm. Pure Appl. Math. 

We match formal asymptotic expansions with differently scaled variables 
to obtain a uniform approximation to the similarity solution of the shock- 
wedge diffraction problem. 


Bokhari, S. H. : Partitioning problems In parallel, pipelined and distributed 

computing. ICASE Report No. 85-54 , November 19, 1985, 34 pages. Submitted 
to IEEE Trans. Comput. 

The problem of optimally assigning the modules of a parallel program 
over the processors of a multiple computer system is addressed. A Sum- 
Bottleneck path algorithm is developed that permits the efficient solution 
of many variants of this problem under some constraints on the structure of 
the partitions. 

In particular, the following problems are solved optimally for a 
single-host, multiple satellite system: Partitioning multiple chain- 

structured parallel programs, multiple arbitrarily structured serial 
programs and single tree structured parallel programs. In addition, the 
problems of partitioning chain structured parallel programs across chain 
connected systems and across shared memory (or shared bus) systems are also 
solved under certain constraints. All solutions for parallel programs are 
equally applicable to pipelined programs. 

These results extend prior research in this area by explicitly 
taking concurrency into account and permit the efficient utilization of 
multiple computer architectures for a wide range of problems of practical 
interest. 


Berger, M. J. and S. H. Bokhari: A partitioning Strategy for non-uniform 

problems on multiprocessors . ICASE Report No. 85-55 , November 20, 1985, 29 
pages. Submitted to IEEE. Trans. Comput. 

We consider the partitioning of a problem on a domain with unequal work 
estimates in different subdomains in a way that balances the work load 
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across multiple processors. Such a problem arises for example in solving 
partial differential equations using an adaptive method that places extra 
grid points in certain subregions of the domain. We use a binary 
decomposition of the domain to partition it into rectangles requiring equal 
computational effort. We then study the communication costs of mapping 
this partitioning onto different multiprocessors: a mesh-connected array, 
a tree machine and a hypercube. The communication cost expressions can be 
used to determine the optimal depth of the above partitioning. 


Bayliss, A., L. Maestrello, P. Parikh, and E. Turkel: Wave, phenomena -in a 

high Reynold* number compile** ible boundary layer. ICASE Report No. 85-56 , 
November 20, 1985, 25 pages. Submitted to Workshop on Stability of Time- 
Dependent and Spatially Varying Flows, (D. L. Dwoyer and M. Y. Hussaini, 
eds.), Springer-Verlag. 

Growth of unstable disturbances in a high Reynolds number compressible 
boundary layer is numerically simulated. Localized periodic surface 
heating and cooling as a means of active control of these disturbances is 
studied. It is shown that compressibility in itself stabilizes the flow 
but at a lower Mach number, significant nonlinear distortions are 
produced. Phase cancellation is shown to be an effective mechanism for 
active boundary layer control. 


Brandt, A. and S. Ta'asan: M ultignld method j $oa nearly *ingular and slightly 
indefinite problem*. ICASE Report No. 85-57 , November 21, 1985, 34 

pages. Lecture Notes in Mathematics, (W. Hackbusch and U. Trottenberg, 
eds.), Springer-Verlag, 1986. 

This paper deals with nearly singular, possibly indefinite problems for 
which the usual multigrid solvers converge very slowly or even diverge. 
The main difficulty is related to some badly approximated smooth functions 
which correspond to eigenfunctions with nearly zero eigenvalues. A 

correction to the usual coarse-grid equations is derived, both in 
Correction Scheme and in Full Approximation Scheme. The performance of the 
new algorithm using this correction is essentially as that of usual 
multigrid for definite problems. 


Burns, J. A. and R. K. Powers: Factorization and reduction method* for 

optimal eoYitnol of distributed parameter *y*tems. ICASE Report No. 85-58 , 
November 22, 1985, 102 pages. Submitted to Mathematics Applicada E 

Computacional . 
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A Chandrasekhar-type factorization method is applied to the linear- 
quadratic optimal control problem for distributed parameter systems. An 
aeroelastic control problem is used as a model example to demonstrate that 
if computationally efficient algorithms, such as those of Chandrasekhar- 
type, are combined with the special structure often available to a 
particular problem, then an abstract approximation theory developed for 
distributed parameter control theory becomes a viable method of solution. 
A numerical scheme based on averaging approximations is applied to 
hereditary control problems. Numerical examples are given. 


Turkel, E., S. Yaniv and U. Landau: Acc.UA.aCi/ ofa *eheme* ^OK. the Entefi 

equation* with non-uni^om me*he*. ICASE Report No. 85-59, December 5, 
1985, 21 pages. AIAA Paper No. 86-0341. 

The effect of non-uniform grids on the solution of the Euler equations 
is analyzed. We consider a Runge-Kutta type scheme based on a finite 
volume formulation. We show that for arbitrary grids the scheme can be 
inconsistent even though it is second-order accurate for uniform grids. An 
improvement is suggested which leads to at least first-order accuracy for 
general grids. Test cases are presented in both two- and three-space 
dimensions. Applications to finite difference and implicit algorithms are 
also given. 


Israeli, M. : UaAching iterative method* j {on. the panabolized and thin tayen. 

Navie^-Stoke* equation*. ICASE Report No. 85-60 , December 5, 1985, 26 
pages. Proc. of Bi-national US/lsrael Workshop on "Impact of 

Supercomputers on the Next Decade of Computational Fluid Dynamics," 
Jerusalem, December 1984. 

Downstream marching iterative schemes for the solution of the 
Parabolized or Thin Layer (PNS or TL) Navier-Stokes equations are 

described. Modifications of the primitive equation global relaxation sweep 
procedure result in efficient second-order marching schemes. These schemes 
take full account of the reduced order of the approximate equations as they 
behave like the SLOR for a single elliptic equation. The improved 
smoothing properties permit the introduction of Multi-Grid acceleration. 
The proposed algorithm is essentially Reynolds number independent and 
therefore can be applied to the solution of the subsonic Euler equations. 
The convergence rates are similar to those obtained by the Multi-grid 

solution of a single elliptic equation; the storage is also comparable as 
only the pressure has to be stored on all levels. Extensions to three- 
dimensional and compressible subsonic flows are discussed. Numerical 

results are presented. 
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Grosch, c. E.: Adapting a Navien.-Stok.e6 code to the. JCL-VAP. ICASE Report 
No. 85-61 , December 1985. Submitted to SIAM J. Sci. Stat. Comput. 

This is a report of the results of an experiment: to adapt a Navier- 
Stokes code, originally developed on a serial computer, to concurrent 
processing on the ICL distributed Array processor (DAP). In this paper the 
algorithm used in solving the Navler-Stokes equations is briefly 
described. The architecture of the DAP and DAP Fortran is also 
described. The modifications of the algorithm so as to fit the DAP are 
given and discussed. Finally, performance results are given and 
conclusions are drawn. 


Canuto, C. and A. Quarteroni: On the. boundany treatment in 6peetn.at method!, 

ion. hypenbotie 6y6tem6. ICASE Report No. 86-1 , December 27, 1985. 

Submitted to J. Comput. Phys. 

Spectral methods have been successfully applied to the simulation of 
slow transients in gas transportation networks. Implicit time advancing 
techniques are naturally suggested by the nature of the problem. 

The aim of this paper is to clarify the correct treatment of the 
boundary conditions in order to avoid any stability restriction originated 
by the boundaries. The Beam and Warming and the Lerat schemes are 
unconditionally linearly stable when used with a Chebyshev pseudospectral 
method. Engineering accuracy for a gas transportation problem is achieved 
at Courant numbers up to 100. 


Gottlieb, D., L. Lustman, and E. Tadmor: Stabiiity anaiy&it, oi 6pe.etn.at 

method:, ion. hypenbotie initiat-boandany vatae 6y&tem . ICASE Report No. 
86-2 , January 16, 1986. Submitted to SIAM J. Numer. Anal. 

We treat a constant coefficient hyperbolic system in one space 
variable, with zero initial data. Dissipative boundary conditions are 
imposed at the two points x = ±1. This problem is discretized by a 
spectral approximation in space. We demonstrate sufficient conditions 
under which the spectral numerical solution is stable — moreover, these 
conditions have to be checked only for scalar equations. The stability 
theorems take the form of explicit bounds for the norm of the solution in 
terms of the boundary data. The dependence of these bounds on N, the 
number of points in the domain (or equivalently the degree of the 
polynomials involved), is investigated for a class of standard spectral 
methods, including Chebyshev and Legendre collocations. 
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Dutt, p.: A R lemann *oloen. ba*ed on a global existence ptioofi h on. the Riemann 
problem. ICASE Report No. 86-3 , January 17, 1986. Submitted to J. Comput. 
Phys . 

Godunov's method and several other methods for computing solutions to 
the equations of gas dynamics use Riemann solvers to resolve 
discontinuities at the interface between cells. A new method is proposed 
here for solving the Riemann problem based on a global existence proof for 
the solution to the Riemann problem. The method is found to be very 
reliable and computationally efficient. 


Saltz, Joel H., Vijay K. Naik, and David M. Nicol: Reduction o£ the ehhect* 

oh the communication delay* In * denticle algorithm* on me** age pa**lng 
MJMV OJidnttectUAe* . ICASE Report No. 86-4 , January 22, 1986, 26 pages. 
Submitted to SIAM J. Sci. Stat. Comput. 

The efficient implementation of algorithms on multiprocessor machines 
requires that the effects of communication delays be minimized. The 
effects of these dealys on the performance of a model problem on a 
hypercube multiprocessor architecture is investigated and methods are 
developed for increasing algorithm efficiency. The model problem under 
investigation is the solution by red-black Successive Over Relaxation 
[Y0UN71] of the heat equation; most of the techniques to be described here 
also apply equally well to the solution of elliptic partial differential 
equations by red-black or multicolor SOR methods. 

This paper identifies methods for reducing communication traffic and 
overhead on a multiprocessor, and reports the results of testing these 
methods on the Intel iPSC Hypercube. We examine methods for partitioning a 
problem's domain across processors, for reducing communication traffic 
during a global convergence check, for reducing the number of global 
convergence checks employed during an iteration, and for concurrently 
iterating on multiple time-steps in a time-dependent problem. Our 
empirical results show that use of these models can markedly reduce a 
numerical problem's execution time. 


Moitra, A., E. Turkel, and A. Kumar: Application oh a Runge-Kutta *ckeme ion 

highspeed inviscid Internal hlout*. ICASE Report No. 86-5 , January 31, 
1986, 26 pages. AIAA Paper NO. 86-0104, AIAA 24th Aerospace Sciences 
Meeting, January 6-9, 1986, Reno, Nevada. 

A multi-stage Runge-Kutta method is analyzed for solving the two- 
dimensional Euler equations for external and internal flow problems. 
Subsonic, supersonic, and highly supersonic flows are studied. Various 
techniques for accelerating the convergence to a steady state are described 
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and analyzed. Effects of the grid aspect ratio on the rate of convergence 
are evaluated. An enthalpy, damping technique applicable to supersonic 
flows is described in detail. Numerical results for supersonic flows 
containing both oblique and normal shocks are presented confirming the 
efficiency of the method. 


Eidson, T. M., M. Y. Hussaini, and T. A. Zang: Simulation o{, the turbulent 

R ayleigh-Benard problem using a s pecXral/ finite deference technique. 
ICASE Report No. 86-6 , February 4, 1986, 43 pages. Proc. EUROMECH 199 
Direct and Large Eddy Simulation of Turbulent Flows, Munich, Germany, 
September 30 - October 2, 1985. 

The three-dimensional, incompressible Navier-Stokes and energy 
equations with the Bousinesq assumption have been directly simulated at a 
Rayleigh number of 3.8 x 10^ and a Prandtl number of 0.76. In the vertical 
direction, wall boundaries were used and in the horizontal, periodic 
boundary conditions were used. A spectral/f inite difference numerical 
method was used to simulate the flow. The flow at these conditions is 
turbulent and a sufficiently fine mesh was used to capture all relevant 
flow scales. The results of the simulation are compared to experimental 
data to justify the conclusion that the small scale motion is adequately 
resolved. 


Nicol, David M. and Paul F. Reynolds, Jr.: An Optimal Repartitioning Decision 

Policy. ICASE Report No. 86-7 , February 10, 1986, 19 pages, IEEE Trans on 
Software Engineering. 

A central problem to parallel processing is the determination of an 
effective partitioning of workload to processors. The effectiveness of any 
given partition is dependent on the stochastic nature of the workload. We 
treat the problem of determining when and if the stochastic behavior of the 
workload has changed enough to warrant the calculation of a new 
partition. We model the problem as a Markov decision process, and derive 
an optimal decision policy. Quantification of this policy is usually 
intractable; we empirically study a heuristic policy which performs nearly 
optimally. Our results suggest that the detection of change is the 
predominant issue in this problem. 


Gottlieb, D., L. Lustman, and E. Tadmor: Convergence o£ Spectral, methods ior 

hyperbolic initial-boundary value systems. ICASE Report No. 86-8 , 
January 17, 1986, 15 pages. Submitted to SIAM J. Numer. Anal. 
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We present here a convergence proof for spectral approximations for 
hyperbolic systems with initial and boundary conditions. We treat in 
detail Chebyshev collocation, but the final result is readily applicable to 
other spectral methods, such as Legendre collocation or tau-methods. 


Abarbanel, S., A. Bayliss, and L. Lustman: Non-reitecting boimdany condition* 

ion. the compressible Navier-Stoke* equation* . ICASE Report No. 86-9 , 
March 4, 1986, 26 pages. Invited lecture at the 7th Internat. Conference 
on Computing Methods in Applied Science and Engineering, December 9-13, 
1985, Versailles, France (INRIA). 

A small perturbation analysis, in the long wavelength regime, is used 
to obtain the downstream boundary condition for the pressure for the flow 
over a flat plate. The methodology is extendable to other geometries. 
Numerical results for high Reynolds number laminar flows show great 
Improvement in convergence rate to steady state as well as in the quality 
of the results. 


Hariharan, S. I.: Nonlinear Acoustic. Wave. Propagation in Atmosphere, ICASE 

Report No. 86-10 , March 5, 1986, 28 pages. IMACS Proceedings on 

Computational Acoustics, August, 1986, Yale. 

In this paper we consider a model problem that simulates an atmospheric 
acoustic wave propagation situation that is nonlinear. The model is 
derived from the basic Euler equations for the atmospheric flow and from 
the regular perturbations for the acoustic part. The nonlinear effects are 
studied by obtaining two successive linear problems in which the second one 
involves the solution of the first problem. Well-posedness of these 
problems is discussed and approximations of the radiation boundary 
conditions that can be used in numerical simulations are presented. 


Banks, H. T. and K. A. Murphy: Quantitative Modeling Growth and Vispersal 
in Population Model* , ICASE Report No, 86-11 , March 1986, 13 pages. 

We discuss techniques for the estimation of nonlinearities and state- 
dependent coefficients in parabolic differential equations. Applications 
to density-dependent population dispersal and nonlinear growth/predation 
models are presented. Computational results using parallel and vector 
architectures are discussed. 
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Lakin, W. D.: A Combined. Intonating and Differentiating Matrix Formulation 

t {on. boundary Value Problem on Rectangular. Domains , icase Report 86-12 , 
March 1986, 30 pages. Submitted to J. Engineering Mathematics. 

Integrating and differentiating matrices allow the numerical 
integration and differentiation of functions whose values are known at 
points of a discrete grid. Previous derivations of these matrices have 
been restricted to one-dimensional grids or to rectangular grids with 
uniform spacing in at least one direction. The present work develops 
integrating and differentiating matrices for grids with non-uniform spacing 
in both directions. The use of these matrices as operators to reformulate 
boundary value problems on rectangular domains as matrix problems for a 
finite-dimensional solution vector is considered. The method requires non- 
uniform grids which include "near-boundary" points. An eigenvalue problem 
for the transverse vibrations of a simply supported rectangular plate is 
solved to illustrate the method. 


Iqbal, M. A., J. H. Saltz, and S. H. Bokhari: Performance tradeoff S in Static 

and dynamic toad balancing strategies. ICASE Report No. 86-13 , March 
1986, 24 pages. Submitted to IEEE Trans. Comput. 

We consider the problem of uniformly distributing the load of a 
parallel program over a multiprocessor system. We analyze a program whose 
structure permits the computation of the optimal static solution. We then 
describe four strategies for load balancing and compare their performance. 

The four strategies are: (1) the optimal static assignment 

algorithm which is guaranteed to yield the best static solution, (2) the 
static binary dissection method which is very fast but sub-optimal, (3) the 
greedy algorithm, a static fully polynomial time approximation scheme, 
which estimates the optimal solution to arbitrary accuracy, and (4) the 
predictive dynamic load balancing heuristic which uses information on the 
precedence relationships within the program and outperforms any of the 
static methods. 

It is also shown that the overhead incurred by the dynamic heuristic 
(4) is reduced considerably if it is started off with a static assignment 
provided by either (1), (2), or (3). 


Turkel , E.: Preconditioned method* for solving the incompressible and tow 

Speed compressibte equations. ICASE No. 86-14 , March 1986, 38 pages. 
Submitted to J. Comp. Phys. 

Acceleration methods are presented for solving the steady state 
incompressible equations. These systems are preconditioned by introducing 
artificial time derivatives which allow for a faster convergence to the 
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steady state. We also consider the compressible equations in conservation 
form with slow flow. Two arbitrary functions ct and g are introduced in 
the general preconditioning. An analysis of this system is presented and 
an optimal value for g is determined given a constant ct • It is further 
shown that the resultant incompressible equations form a symmetric 
hyperbolic system and so are well-posed. Several generalizations to the 
compressible equations are presented which .generalize previous results. 


Gibson, J. S. and I. G. Rosen: Hwme>vieaZ approximation faoh. the infinite- 

dimensional du> arete-time optimal tinear-qaadratia regulator problem. 
ICASE Report No. 86-15 , March 20, 1986, 82 pages. Submitted to SIAM J. 
Cont. Opt. 

An abstract approximation framework is developed for the finite and 
infinite time horizon discrete-time linear-quadratic regulator problem for 
systems whose state dynamics are described by a linear semigroup of 
operators of an infinite dimensional Hilbert space. The schemes included 
in the framework yield finite dimensional approximations to the linear 
state feedback gains which determine the optimal control law. Convergence 
arguments are given. Examples involving hereditary and parabolic systems 
and the vibration of a flexible beam are considered. Spline-based finite 
element schemes for these classes of problems, together with numerical 
results, are presented and discussed. 


Gibson, J. S. and I. G. Rosen: Shifting -the closed-loop spectrum in the 

optimal linear quadratic regulator problem for hereditary systems . icase 
R eport No. 86-16 , March 20, 1986, 47 pages. Submitted to IEEE Trans. Auto 
J. Control. 

In the optimal linear quadratic regulator problem for finite 
dimensional systems, the method known as an a-shift can be used to 
produce a closed-loop system whose spectrum lies to the left of some 
specified vertical line; that is, a closed-loop system with a prescribed 
degree of stability. This paper treats the extension of the a-shift to 
hereditary systems. As in finite dimensions, the shift can be accomplished 
by adding a times the identity to the open-loop semigroup generator and 
then solving an optimal regulator problem. However, this approach does not 
work with a new approximation scheme for hereditary control problems 
recently developed by Kappel and Salamon. Since this scheme is among the 
best to date for the numerical solution of the linear regulator problem for 
hereditary systems, an alternative method for shifting the closed-loop 
spectrum is needed. An a-shift technique that can be used with the 
Kappel-Salamon approximation scheme is developed. Both the continuous-time 
and discrete-time problems are considered. A numerical example which 
demonstrates the feasibility of the method is included. 
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Shih, Ke-Gang: On the existence o & 4 olutions &or an aquation arising in the. 

theory laminar &low in a uniformly porous channel with injection, icase 
R eport No. 86-17 , March 21, 1986, 26 pages. Submitted to SIAM J. Appl. 
Math. 


We establish the existence of concave solutions of Berman's equation 
which describes the laminar flow in channels with injection through porous 
walls. It is found that the (unique) concave solutions exist for all 
injection Reynolds number R < 0. 


South, Jr., J. C. and M. Y. Hussaini, (eds.): Advance* in Numerical and 

Applied Mathematics, ICASE Report No. 86-18 , March 28, 1986, to appear in 
Applied Numerical Mathematics. 

This collection of papers covers some recent developments in numerical 
analysis and computational fluid dynamics. Some of these studies are of a 
fundamental nature. They address basic issues such as intermediate 
boundary conditions for approximate factorization schemes, existence and 
uniqueness of steady states for time-dependent problems, pitfalls of 
implicit time stepping, etc. The other studies deal with modern numerical 
methods such as total variation-diminishing schemes, higher-order variants 
of vortex and particle methods, spectral multidomain techniques, and front- 
tracking techniques. There is also a paper on adaptive grids. The fluid 
dynamics papers treat the classical problems of incompressible flows in 
helically-coiled pipes, vortex breakdown, and transonic flows. 


Lee, Tai-yong, Costas Kravaris, and John H. Seinfeld; History Matching by 
Spline Approximation and Regularization in Single-Phase Areal Reservoirs. 
ICASE Report No. 86-19 , March 31, 1986, 53 pages. Submitted to SPE 

Reservoir Engineering. 

An automatic history matching algorithm is developed based on bi-cubic 
spline approximations of permeability and porosity distributions and on the 
theory of regularization to estimate permeability or porosity in a single- 
phase, two-dimensional areal reservoir from well pressure data. The 
regularization feature of the algorithm, the theoretical details of which 
are described by Kravaris and Seinfeld is used to convert the ill-posed 
history matching problem into a well-posed problem. The algorithm employs 
the conjugate gradient method of Nazareth as its core minimization 
method. A number of numerical experiments are carried out to evaluate the 
performance of the algorithm. Comparisons with conventional (non- 
regularized) automatic history matching algorithms indicate the superiority 
of the new algorithm with respect to the parameter estimates obtained. A 
quasioptimal regularization parameter is determined without requiring a 
priori information on the statistical properties of the observations. 



ICASE COLLOQUIA 


October 


October 


October 


November 


November 


November 


December 


December 


December 


December 


December 


1 Hr. Gary Schmidt, Mr. Henno Allik and Mr. Bob Thomas, Bolt, 

Beranek and Newman, Inc.: Ove> ivlew and Use of the Butterfly 

System 


7 Professor James Glimm, Courant Institute of Mathematical 

Sciences: An Overview o f the Front Tracking Method 


24 Professor Victor Takhot, Princeton University: 

Renormallzatlon-Group Based Large Eddy and Turbulence 
Tram port Models 


1 Professor David Casasent, Carnegie-Mellon University: 

Optical Computing 


20 Professor Roddam Narasimha, California Institute of 

Technology: The Transition Zone In the Boundary Layer 


27 Dr. Sanjay M. Correa, General Electric Corporate Research 

and Development Center: Models for Turbulent Diffusion 

Flames 


3 Professor Youcef Saad, Yale University: Linear Algebra 

Algorithms on the Hypercube 

5 Professor Mel S. Berger, M.I.T.: Vortices In Three- 

Vlmemlonal Ideal Fluids as a Sollton Theory 


11 Professor David Kopriva, Florida State University: A 

Multidomain Spectral Method for Hyperbolic Equatlom 


12 Professor Mel S. Berger, M.I.T.: Vortices In Modern Physics 

(An Example) 


13 Professor Lawrence Sirovich, Brown University: Neu) 

Eigenfunction Methods with Appllcatlom to Turbulence and 
Other Pattern Recognition Problems 


42 



December 


December 


January 


January 


January 


January 


February 


February 


February 


February 


February 


43 


19 Dr. Alvin Bayliss, Exxon Corporate Research: Numerical 

Analysis Ojj a Hoph Bihurcation Using an Adaptive. Pseudo- 
Spectral Method 

20 Professor Mel S. Berger, M.I.T.: Bihurcation and Complete 

Integrability h ofL hi online an. Dynamical Systems 


13 Professor Alton Highsmith, Texas A & M University: The 

Damage Induced Stress Redistribution in Composite Laminates 


16 Mr. James B. Salem, Thinking Machines Corporation: The 

Connection Machine : Massively Parallel Computing and 

Modeling Fluid Flow with Cellulan Automata 

21 Professor Ivo Babuska, University of Maryland: The h- p 

tension oh the Finite Element Method. Theory and Practice 


29 Professor J. L. Lions, President, CNES: Boundary 

Stabilization and Controllability oh Hyperbolic and 
Petrowsky Systems 


4 Professor Paul F. Reynolds, Jr., University of Virginia: 

Static Incremental Detection oh Synchronization Errors in 
Parallel Programs 


6 Dr. Thomas D. Taylor, Dr. Hwar C. Ku, The John Hopkins 

Applied Physics Laboratory: Pseudospectral Methods and 

Boundary Conditions h ofL Solution oh the Incompressible 
Navier-Stokes Equations 


18 Dr. Eitan Tadmor, Tel-Aviv University and ICASE: The 

Numerical Viscosity oh Entropy Stable Schemes h 0>l Systems oh 
Conservation Laws 


20 Mr. David Kamowitz, University of Wisconsin - Madison: 

Theoretical and Computational Results h ofL MGR Multigrid 
Methods 


25 Mr. Stephen L. Keeling, University of Tennessee: 

Galerkin/Runge-Kutta Discretizations h on - Parabolic Partial 
Dihh&rential Equations 



44 


March 


March 


March 


4 Mr. Anthony T. Chronopoulis , University of Illinois: 

Conjugate Gradient- Like. Algorithm ion. Solving Nonline an. 
Systems o& Equations Arising in the Numenieal Integration oi 
StiU OVE’s 


19 Dr. William J. Feiereisen, Brown Boveri & Company: The 

Modelling o& Transition and Surface Roughness Ejects in 
Boundary -Layer Flows 


25 Professor George Adomian, The University of Georgia: 

Analytical Solution oi Nonlinear Systems 



ICASE STAFF 


I. ADMINISTRATIVE 

Milton E. Rose, Director (Through December 1985) 

Ph.D., Mathematics, New York University, 1953 
Numerical Methods 

Robert G. Voigt, Director (Beginning January 1986) 

Ph.D., Mathematics, University of Maryland, 1969 
Numerical and Computational Techniques 

Linda T. Johnson, Administrative Assistant 

Etta M. Blair, Personnel/Bookkeeping Secretary 

Barbara A. Cardasis, Administrative Secretary 

Sidney A. Chappell, Technical Publications/Summer Housing Secretary 
(Beginning March 24, 1986) 

Carla J. Hult, Office Assistant (Beginning February 21, 1986) 

Barbara A. Kraft, Senior Technical Publications Secretary 

Barbara A. Rohrbach, Technical Publications Secretary (Beginning January 22, 
1986) 

Emily N. Todd, Visitor Coordinator/Correspondence Secretary 

Georgia A. Voss, Technical Publications/Summer Housing Secretary (Through 
February 21, 1986) 


II. SCIENCE COUNCIL for APPLIED MATHEMATICS and COMPUTER SCIENCE 


Bruce Arden, Chairman and Arthur Doty, Professor, Department of Electrical 
Engineering and Computer Science, Princeton University. 

Andrew J. Callegari, Director, Theoretical & Mathematical Sciences Laboratory, 
Exxon Research & Engineering Company. 

Peter Denning, Director, RIACS, NASA/ Ames Research Center. 

Michael J. Flynn, Professor, Department of Electrical Engineering, Computer 
Systems Laboratory, Stanford University. 

Bernard Galler, Professor, Department of Computer and Communication Sciences 
and Associate Director of the Computer Center, University of Michigan. 


45 



46 


Anthony C. Hearn, Department Head, Department of Information Sciences, Rand 
Corporation. 

Herbert Keller, Professor, Physics, Math, and Astronomy, California Institute 
of Technology. 

Seymour V. Parter, Professor, Department of Mathematics, University of 
Wisconsin. 

Werner C. Rheinboldt, Andrew W. Mellon Professor, Department of Mathematics 
and Statistics, University of Pittsburgh. 

Milton E. Rose, Director, Institute for Computer Applications in Science and 
Engineering, NASA Langley Research Center. (Through December 1985) 

Robert G. Voigt, Director, Institute for Computer Applications in Science and 
Engineering, NASA Langley Research Center. (Beginning January 1986) 


III. ASSOCIATE MEMBERS 


Saul S. Abarbanel, Professor, Department of Applied Mathematics, Tel-Aviv 
University, Israel. 

H. Thomas Banks, Professor, Division of Applied Mathematics, Brown University. 

Peter D. Lax, Professor, Courant Institute of Mathematical Sciences, New York 
University. 

William R. Sears, Professor, Department of Aerospace and Mechanical 
Engineering, University of Arizona. 


IV. SENIOR STAFF SCIENTISTS 


David Gottlieb - Ph.D., Numerical Analysis, Tel-Aviv University, Israel, 
1972. Professor, Department of Applied Mathematics, Brown University. 
Numerical Methods for Partial Differential Equations. (July 1974 to January 
1987) 

M. Yousuff Hussaini - Ph.D., Mechanical Engineering, University of California, 
1970. Computational Fluid Dynamics. (Beginning April 1978) 



47 


V. SCIENTIFIC STAFF 


Pravir Dutt - Ph.D., Mathematics, University of California, Los Angeles, 
1985. Numerical Analysis related to Navier-Stokes Equations. (July 1985 to 
June 1987) 

Thomas Jackson - Ph.D., Mathematics, Rensselaer Polytechnic Institute, 1985. 
Numerical and Analytical Methods for Chemically Reacting Flows. (October 1985 
to October 1987) 

Liviu Lustman - Ph.D., Applied Mathematics, Tel-Aviv University, Israel, 
1979. Numerical Methods for Partial Differential Equations, Spectral Methods, 
Computational Fluid Dynamics. (August 1984 to May 1986) 

Nessan Mac Giolla Mhuiris - Ph.D., Theoretical and Applied Mechanics, Cornell 
University, 1986. Theoretical and Numerical Models of Vortex Breakdown. 
(October 1985 to October 1987) 

Vi jay Naik - M.S., Mechanical Engineering, University of Florida, 1982. AM., 
Computer Science, Duke University, 1984. Parallel Processing, Computational 
Fluid Dynamics. (August 1984 to August 1986) 

David M. Nicol - Ph.D., Computer Science, University of Virginia, 1985. 
Techniques for Mapping Algorithms onto Parallel Computing Systems. (September 
1985 to August 1987) 

Douglas Peterson - B.S., Computer Science, Christopher Newport College, 
1984. Computer Systems. (October 1985 to October 1987) 

Robert K. Powers - Ph.D., Applied Mathematics, Virginia Polytechnic Institute 
and State University, 1984. Numerical Techniques for Optimal Control. (March 

1984 to September 1986) 

Georg Propst - Ph.D., Mathematics-Functional Analysis, Karl Franzens 
Universitat, Austria, 1985. Linear Quadratic Control Problem for Delay 
Systems. (September 1985 to August 1987) 

Joel H. Saltz - Ph.D., Computer Science, Duke University, 1985. Analysis of 
Parallel Methods for Time Dependent Partial Differential Equations. (July 

1985 to September 1986) 

Ke-Gang Shih - Ph.D., Applied Mathematics, University of Maryland, 1985. 
Analysis of Methods to Treat Problems in Nonlinear Elasticity. (August 1985 
to July 1987) 

Nancy E. Shoemaker - Ph.D., Mathematics, Rensselaer Polytechnic Institute, 
1977. Computer Systems, particularly UNIX. (February 1982 to January 1986) 

Shlomo Ta'asan - Ph.D., Applied Mathematics, Weizman Institute of Science, 
Israel, 1984. Multigrid Methods. (October 1984 to September 1987) 



48 


VI. VISITING SCIENTISTS 


Saul S. Abarbanel - Ph.D., Theoretical Aerodynamics, Massachusetts Institute 
of Technology, 1959. Professor, Department of Applied Mathematics, Tel-Aviv 
University, Israel. Numerical Analysis of Partial Differential Equations. 
(January to December 1986) 

David C. Arney - Ph.D., Mathematics, Rensselaer Polytechnic Institute, 1985. 
Associate Professor, Department of Mathematics, United States Military 
Academy, West Point. (September 1985 to January 1986) 

Melvyn S. Berger - Ph.D., Mathematics, Yale University, 1964. Visiting 
Professor, Department of Mathematics, Massachusetts Institute of Technology. 
Partial Differential Equations. (December 1985 to January 1986) 

Shahid H. Bokhari - Ph.D., Electrical and Computer Engineering, University of 
Massachusetts, Amherst, 1978. Associate Professor, Department of Electrical 
Engineering, University of Engineering and Technology, Lahore, Pakistan. 
Parallel Processing, Distributed Computing and Computer Architecture. (July 
1984 - June 1986) 

Eitan Tadmor - Ph.D., Numerical Analysis, Tel-Aviv University, 1979. Senior 
Lecturer, Department of Applied Mathematics, Tel-Aviv University, Israel. 
Numerical Methods for Partial Differencial Equations. (February 1986) 

Hillel Tal-Ezer - Ph.D., Applied Mathematics, Tel-Aviv University, 1985. 
Ph. D. , Department of Mathematics, Tel-Aviv University, Israel. Spectral 
Methods for Partial Differential Equations. (February 1986) 

Eli Turkel - Ph.D., Applied Mathematics, New York University, 1970. Associate 
Professor, Department of Applied Mathematics, Tel-Aviv University, Israel. 
(January - December 1985) 


VII. CONSULTANTS 


Loyce M. Adams - Ph.D., Applied Mathematics, University of Virginia, 1983. 
Parallel Numerical Algorithms. Associate Professor, Department of Applied 
Mathematics, University of Washington. Numerical Methods for Parallel 
Computing Systems. 

Ivo Babuska - Ph.D., Technical Sciences, Tech University, Prague, 
Czechoslovakia, 1952; Mathematics, Academy of Sciences, Prague, 1956; D.Sc., 
Mathematics, Academy of Science, Prague, 1960. Professor, Institute for 
Physical Science and Technology, University of Maryland. Numerical Methods 
for Partial Differential Equations. 



49 


Alvin Bayliss - Ph.D., Mathematics, New York University, 1975. Senior Staff 
Scientist, Exxon Corporate Research. Numerical Methods for Partial 
Differential Equations. 

Marsha J. Berger - Ph.D., Numerical Analysis, Stanford University, 1982. 
Research Associate, Courant Institute of Mathematical Sciences. Numerical 
Methods for Partial Differential Equations. 

Dennis W. Brewer - Ph.D., Mathematics, University of Wisconsin, Madison, 
1975. Associate Professor, Department of Mathematical Sciences, University of 
Arkansas. Methods for Parameter Identification and Estimation. 

John A. Burns - Ph.D., Mathematics, University of Oklahoma, 1973. Professor, 
Virginia Polytechnic Institute and State University. Numerical Methods in 
Feedback Control and Parameter Estimation. 

James M. Crowley - Ph.D., Applied Mathematics, Brown University, 1982. 
Assistant Professor, Department of Mathematics, (Captain), U.S. Air Force 
Academy. Identification for Distributed Parameter Systems. 

Janice E. Cuny - Ph.D., Computer Science, University of Michigan, 1981. 
Assistant Professor, Department of Computer and Information Science, 
University of Massachusetts. Program Environments for Parallel Computing 
Systems . 

Peter R. Eiseman - Ph.D., Mathematics, University of Illinois, 1970. Senior 
Research Scientist and Adjunct Professor, Department of Applied Physics and of 
Nuclear Engineering, Columbia University. Computational Fluid Dynamics. 

Stefan Feyock - Ph.D., Computer Science, University of Wisconsin, 1971. 
Associate Professor, Department of Mathematics and Computer Science, College 
of William and Mary. Artificial Intelligence. 

George J. Fix - Ph.D., Mathematics, Harvard University, 1968. Professor and 
Chairman, Department of Mathematics, Carnegie-Mellon University. Numerical 
Methods for Partial Differential Equations. 

Joseph E. Flaherty - Ph.D., Applied Mechanics, Polytechnic Institute of 
Brooklyn, 1969. Professor and Chairman, Departments of Mathematical Sciences 
and Computer Science, Rensselaer Polytechnic Institute. Singular 
Perturbations, Numerical Analysis, and Applied Mathematics. 

Dennis B. Gannon - Ph.D., Computer Science, University of Illinois, 1980. 
Assistant Professor, Department of Computer Science, Indiana University at 
Bloomington. Numerical Methods, Software and Computer Architecture Design. 

James F. Geer - Ph.D., Applied Mathematics, New York University, 1967. 
Professor, Systems Science and Mathematical Sciences, Watson School of 
Engineering, Applied Science and Technology, SUNY-Binghamton. Perturbation 
Methods and Asymptotic Expansions of Solutions to Partial Differential 
Equations. 



50 


J. Steven Gibson - Ph.D., Engineering Mechanics, University of Texas at 
Austin, 1975. Associate Professor, Department of Mechanical, Aerospace and 
Nuclear Engineering, University of California at Los Angeles. Control of 
Distributed Systems. 

Maya Gokhale - Ph.D., Computer Science, University of Pennsylvania, 1983. 
Assistant Professor, Department of Computer Science, University of Delaware. 
Nonprocedural Languages and Parallel Processing. 

Chester E. Grosch - Ph.D., Physics - Fluid Dynamics, Stevens Institute of 
Technology, 1967. Professor, Department of Computer Science and Slover 
Professor, Department of Oceanography, Old Dominion University. Hydrodynamic 
Stability, Computational Fluid Dynamics, Unsteady Boundary Layers and 
Algorithms for Array Processors. 

Max D. Gunzburger - Ph.D., Mathematics, New York University, 1969. Professor, 
Department of Mathematics, Carnegie-Mellon University. Fluid Dynamics and 
Numerical Methods for Partial Differential Equations. 

Subramaniya I. Hariharan - Ph.D., Mathematics, Carnegie-Mellon University, 

1980. Associate Professor, Department of Mathematical Sciences, University of 
Akron. Numerical Methods for Partial Differential Equations. 

Amiram Harten - Ph.D., Mathematics, New York University, 1974. Associate 
Professor, Department of Mathematics, Tel-Aviv University, Israel. Numerical 
Solution for Partial Differential Equations. 

Kazufumi Ito - Ph.D., Systems Science and Mathematics, Washington University, 

1981. Assistant Professor, Division of Applied Mathematics, Brown 
University. Control Theory. 

Harry F. Jordan - Ph.D., Physics, University of Illinois, 1977. Professor 
Department of Electrical and Computer Engineering, University of Colorado at 
Boulder. Parallel Computation. 

Ashwani K. Kapila - Ph.D., Theoretical and Applied Mechanics, Cornell 
University, 1975. Associate Professor, Department of Mathematical Sciences, 
Rensselaer Polytechnic Institute. Ordinary and Partial Differential 
Equations, Asymptotic Methods. 

David A. Kopriva - Ph.D., Applied Mathematics, University of Arizona, 1982. 
Assistant Professor, Department of Mathematics, Florida State University. 
Spectral Methods for Problems in Fluid Dynamics. 

Heinz-Otto Kreiss - Ph.D., Mathematics, Royal Institute of Technology, Sweden 
1960. Professor, Department of Applied Mathematics, California Institute of 
Technology. Numerical Analysis. 

Karl K. Kunish - Ph.D., Mathematics, University of Graz, Austria. Associate 
Professor, Department of Mathematics, Technical University of Graz, Austria. 
Parameter Identification and Control. 



51 


William D. Lakin - Ph.D., Applied Mathematics, University of Chicago, 1968. 
Eminent Professor, Department, of Mathematical Sciences, Old Dominion 
University. Fluid Mechanics and Elastic Vibrations. 

Patricia K. Lamm - Ph.D., Applied Mathematics, Brown University, 1981. 
Assistant Professor, Department of Mathematics, Southern Methodist University. 
Control and Identification of Partial Differential Equations. 

Robert W. MacCormack - M.S., Mathematics, Stanford University. Professor, 
Department of Aeronautics and Astronautics, Stanford University. 
Computational Fluid Dynamics and Numerical Analysis. 

Piyush Mehrotra - Ph.D., Computer Science, University of Virginia, 1982. 
Associate Professor, Department of Computer Science, Purdue University. 
Programming Languages for Multiprocessor Systems. 

Katherine A. Murphy - Ph.D., Applied Mathematics, Brown University, 1983. 
Visiting Assistant Professor, Department of Applied Mathematics, Brown 
University. Control Theory and Estimation of Parameters. 

Roy A. Nicolaides - Ph.D., Mathematics, University of London, 1972. 
Professor, Department of Mathematics, Carnegie-Mellon University. Numerical 
Solution of Partial Differential Equations. 

Robert Noonan - Ph.D., Computer Science, Purdue University, 1971. Professor, 
Department of Computer Science, College of William and Mary. Parallel 
Programming Techniques. 

Joseph E. Oliger - Ph.D., Numerical Analysis, Uppsala University, Sweden, 
1973. Associate Professor, Department of Computer Science, Stanford 
University. Numerical Methods for Partial Differential Equations. 

James E. Ortega - Ph.D., Mathematics, Stanford University, 1962. Professor 
and Chairman, Department of Applied Mathematics, University of Virginia. 
Numerical Methods for Partial Differential Equations. 

Stanley Osher - Ph.D., Functional Analysis, New York University, 1966. 
Professor, Department of Mathematics, University of California at Los Angeles. 
Methods for the Numerical Analysis of Partial Differential Equations. 

Anthony Patera - Ph.D., Applied Mathematics, Massachusetts Institute of 
Technology, 1982. Assistant Professor, Department of Mechanical Engineering, 
Massachusetts Institute of Technology. Numerical Methods for Fluid Dynamics 
and Heat Transfer. 

Merrell L. Patrick - Ph.D., Mathematics, Carnegie-Mellon University, 1964. 
Professor, Department of Computer Science, Duke University. Parallel 
Numerical Algorithms and Architectures. 

Terrence W. Pratt - Ph.D., Mathematics/Computer Science, University of Texas 
at Austin, 1965. Professor, Department of Computer Science, University of 
Virginia. Programming Languages. 



52 


Daniel A. Reed - Ph.D., Computer Science, Purdue University, 1983. Assistant 
Professor, Department of Computer Science, University of Illinois. Parallel 
Processing. 

Paul F. Reynolds - Ph.D., Computer Science, University of Texas at Austin, 
1979. Assistant Professor, Department of Computer Science, The University of 
Virginia. Parallel Computing Systems. 

I. Gary Rosen - Ph.D., Applied Mathematics, Brown University, 1980. Assistant 
Professor, Department of Applied Mathematics, University of Southern 
California. Numerical Methods for Problems in Control Systems. 

Nancy E. Shoemaker - Ph.D., Mathematics, Rensselaer Polytechnic Institute, 
1977. Utilization of the Computational Resources. 

Lawrence Sirovich - Ph.D., Fluid Mechanics, Johns Hopkins University, I960.. 
Professor, Division of Applied Mathematics, Brown University. Fluid Mechanics. 

Charles G. Speziale - Ph.D., Aerospace and Mechanical Sciences, Princeton 
University, 1978. Associate Professor, School of Mechanical Engineering, 
Georgia Institute of Technology. Turbulence Modeling. 

Edward A. Spiegel - Ph.D, Astronomy, University of Michigan, 1958. Rutherford 
Professor of Astronomy, Department of Astronomy and Astrophysics, Columbia 
University. Mathematical Aspects of Turbulence. 

John C. Strikwerda - Ph.D., Mathematics, Stanford University, 1976. Associate 
Professor, Department of Computer Sciences, University of Wisconsin-Madison. 
Numerical Methods for Partial Differential Equations. 

J. Christian Wild - Ph.D., Computer Science, Rutgers University, 1977. 
Assistant Professor, Department of Computer Science, Old Dominion 
University. Concurrent Computing Systems. 

John Van Rosendale - Ph.D., Computer Science, University of Illinois, 1980. 
Assistant Professor, Department of Computer Science, University of Utah. 
Parallel Processing. 


VIII. STUDENT ASSISTANTS 


Mary Ann 0. Bynum - Graduate student at the College of William and Mary. 
(September 1984 to Present) 

Raad A. Fatoohi - Graduate student at Old Dominion University. (September 
1985 to Present) 



53 


IX. GRADUATE FELLOWS 

Ashraf M. Iqbal - Graduate student at University of Engineering and 
Technology, Lahore, Pakistan. (April 1985 to Present) 



Standard Bibliographic Page 


1. Report No. 2. Government Accession No. 

NASA CR- 178 119 

3. Recipient’s Catalog No. 

4. Title and Subtitle 

Semi-Annual Report 

October 3, 1985 through March 31, 1986 

5. Report Date 

April 1986 

6. Performing Organization Code 

7. Author(s) 

8. Performing Organization Report No. 

9. Performing Organization Name and Address 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23665 

10. Work Unit No. 

11. Contract or Grant No. 
NAS1-18107 

13. Type of Report and Period Covered 
Contractor Report 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 

505-31-83-01 


15. Supplementary Notes 


Langle Technical Monitor: J. C. South 


Final Report 


16. Abstract 


This report summarizes research conducted at the Institute for Computer 
Applications in Science and Engineering in applied mathematics, numerical 
analysis, and computer science during the period October 3, 1985 through 
March 31, 1986. 


17. Key Words (Suggested by Authors(s)) 

Mathematics 
Numerical Analysis 
Computer Science 


18. Distribution Statement 


59 - Mathematics and Computer Sciences 
(General) 

Unclassified - Unlimited 


19. Security Classif.(of this report) 

20. Security Classif.(of this page) 

21. No. of Pages 

22. Price 

Unclassified 

Unclassified 

57 

A0 4 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1986 




















End of Document 



